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Abstract 

We consider the possibility of neutrinos being Dirac particles and study minimal mass 
matrices with as much zero entries as possible. We find that up to 5 zero entries are 
allowed. Those matrices predict one vanishing mass state, CP conservation and Ues 
either zero or proportional to R, where R is the ratio of the solar and atmospheric 
Am^. Matrices containing 4 zeros can be classihed in categories predicting Ue 3 = 0, 
Ue 3 / 0 but no CP violation or |? 7 e 3 | 7 ^ 0 and possible CP violation. Some cases allow 
to set constraints on the neutrino masses. The characteristic value of Ues capable 
of distinguishing some of the cases with non-trivial phenomenological consequences 
is about R/2 sin 20 i 2 . Matrices containing 3 and less zero entries imply (with a 
few exceptions) no correlation for the observables. We outline models leading to 
the textures based on the Froggatt-Nielsen mechanism or the non-Abelian discrete 
symmetry D 4 x Z 2 . 
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1 Introduction 


Recent years witnessed a dramatic increase of our knowledge about the properties of neu¬ 
trinos (for recent reviews, see for instance P). The experimental evidence showed that at 
least two neutrinos possess a non-vanishing rest mass, whose value is bounded from above 
by roughly 1 eV This implied smallness of neutrino masses is commonly attributed 

to the see-saw mechanism jl], which explains the tiny neutrino masses through their in¬ 
verse proportionality to the mass scale of heavy Majorana neutrinos. In addition, the 
light neutrinos are predicted to be of Majorana type. The smoking-gun signature of this 
property would be the observation of neutrinoless double beta decay [3j. There is currently 
no incontrovertible evidence for this process (see the discussion in (3]) and therefore the 
possibility that neutrinos are Dirac particles, just like all the other known fermions, has to 
be considered as well. If neutrinos are indeed Dirac particles, their mass term is given just 
by the coupling of the lepton doublets with the Higgs doublet and a right-handed neutrino 
Standard Model singlet state. One should then explain however why the corresponding 
Yukawa coupling is so much smaller than the Yukawa coupling for the other fermions. 
There are several possibilities for this to arise. All have in common that effectively there 
is an additional D(1)b-l symmetry which forbids a Majorana mass term. For instance: 

• in theories with large extra dimensions [5] the effective Yukawa coupling is determined 
by the overlap of the wave functions of the SM particles (residing in the usual 3 -|- 
1 dimensions) with the wave functions of the right-handed singlet neutrinos. The 
latter propagate in additional dimensions (the bulk) and a volume suppression of 
order Mp/M-pi occurs, where Mpi is the Planck scale and Mp the fundamental gravity 
scale, possibly as low as 10 TeV. Similar volume suppression can occur in models 
with warped extra dimensions or “fat brane” models [B]. For more details of these 
possibilities, see the recent reviews in [7]; 

• in supersymmetric frameworks the suppression factor might be given by the SUSY 
breaking scale divided by the Planck scale in scenarios with anomaly mediated su¬ 
persymmetry breaking jH] . For other supersymmetry or supergravity models of Dirac 
neutrino masses, see, e.g., jU]; 

• in superstring models it is possible to obtain small neutrino masses through a hierar¬ 
chy of vacuum expectation values of SM singlets UDI. In this respect it is interesting 
to note that in string theories it seems difficult to obtain a simple see-saw mechanism 

ffH; 

• if neutrinos are massless at tree level, small Dirac masses might be generated through 
radiative mechanisms, so that the smallness is explained by loop suppression [T^ : 

• other approaches to generate Dirac neutrino masses are discussed in Refs. ng. 

A common argument in favor of Majorana neutrinos is that the implied lepton number vi¬ 
olation can be used to generate the baryon asymmetry of the universe via the leptogenesis 
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mechanism However, even with Dirac neutrinos it is possible to achieve this usmn]. 
We therefore conclude that there is some motivation for studying Dirac neutrinos. 

Once the Dirac nature of neutrinos is accepted, one can ask immediately what the struc¬ 
ture of the neutrino mass matrix rriy looks like. In particular, minimal mass matrices are 
of interest since they usually imply interesting and testable correlations between the ob¬ 
servables. For Majorana neutrinos, there has been some interest in texture zeros of 
HU- It has been found that only two independent zero entries are allowed by current data 
HU- Models leading to the successful textures HHI as well as renormalization effects HSI 
have been studied. Here we perform a similar analysis for Dirac neutrino mass matrices. 
Since Dirac mass matrices are — in contrast to Majorana mass matrices — in general not 
symmetric, the results are expected to differ. We shall hnd that one is able to put up to 5 
zero entries in a Dirac neutrino mass matrix. Those cases predict one vanishing neutrino 
mass, CP conservation and small or zero f/gs. Already for matrices with 4 zeros there are 
successful cases in which there is no correlation between the observables. Some of them 
display however a correlation different from the 5 zero cases. 

We also outline models based on the discrete symmetry group D 4 x Z 2 or on the Froggatt- 
Nielsen mechanism, which allow in principle to reproduce the textures under consideration. 

The paper is structured as follows: Section El contains several general considerations about 
Dirac mass matrices, their form as implied by current data and comment on the possibility 
to determine the nature of the neutrino. Then we perform a systematic search for zero 
textures for Dirac neutrinos in Sectional discussing the cases of 5, 4 and 3 or less zeros. We 
comment on the implications of the Dirac mass matrices being symmetric. Models based 
on the flavor symmetry D 4 x Z 2 and on the Froggatt-Nielsen mechanism are outlined in 
Section m before we conclude in Section 0 


2 Dirac Neutrino Mass Matrices: General Consider¬ 
ations 


For Dirac neutrinos, the PMNS (20] niatrix can be parametrized as 



/ 

C12C13 

■S12C13 

Si3 

u = 


-S12C23 - C12S23S13 e“*'^ 

C12C23 - S12S23S13 

S23C13 


V 

S12S23 - C12C23S13 

-C12S23 - S12C23S13 e“*'^ 

C23C13 


( 1 ) 


where Cij = cosOij, Sij = sin 6 'ij and 5 denotes the CP phase, sometimes called “Dirac- 
phase”. Global analyzes of current data allow for the following values of the observables 
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sin^ 012 = 0.30 +o;o5 , 


sin^ 023 = 0.50 +o;i2 , 


sin^0i3 = |t43p < 0.031 

Am^i = AttIq = (T.O^Q^g) ■ 10“® eV^ , = Am^ = (2.2 ^q; 5) ■ 10“^ eV^ ^ 2 ^ 

p ^ ^^0 _ Q 036 +°-°^'^ 

^ “ A^ “ U.03b_o.oii , 

where Am^ denotes rrif — m|. We gave the best-ht values and the allowed 2 ct ranges. 
Regarding the mass ordering, both the normal and the inverted ordering can be accommo¬ 
dated, corresponding to the cases m 3 > m 2 > mi and m 2 > mi > m 3 , respectively. The 
extreme cases are the normal hierarchy (NH), which is dehned as m| ~ Am\ ^ AmQ ~ 
m^ 3> mf, and the inverted hierarchy (IH), which is dehned as m^ ~ m^ ~ Am\ ^ m|. 
The quasi-degenerate spectrum (QD) occurs when the common neutrino mass scale mo is 
much larger than the mass splittings, i.e., m^ = m^ ~ m^ — ml ^ Am^ 


Mass matrices for Dirac neutrinos are not symmetric. Hence, they are parameterized by 
nine complex variables a,... ,1: 

( a b d \ 

e f g ] =17mfsRt ^ ( 3 ) 

h k I ) 


where f/, V are unitary matrices and is a diagonal mass matrix containing the three 
mass eigenstates. In the basis in which the charged lepton mass matrix is diagonal the 
relevant Lagrangian reads 

C = iR + -^W^iL'y^{i^i)L + h.c. , (4) 

where i = e,/i, r. After substituting = U and = V (i'e)R we identify the 
PMNS matrix U. This matrix is associated with the diagonalization of m^m^, i.e., 

U'^mym^U = hU = . (5) 

Here we also dehned the matrix h = m^ml. The matrix V diagonalizes mlm^ and is, 
in the absence of right-handed currents, unobservable. In case of Majorana neutrinos, we 
would have V = U*. 


One of the most interesting aspects of neutrino physics is the possibility of leptonic CP 
violation. Whether a given neutrino mass matrix implies CP violation in oscillation ex¬ 
periments (note that due to the Dirac nature assumed in this work there are no so-called 
Majorana phases j221) can for instance be determined from the characteristic matrix h. In 
terms of h, one can write the Jarlskog invariant Jcp as m 


J = Im(hi 2 0-23 0 - 31 ) 
Am^i Am^.^ Arn^2 


Am 

i sin 2012 sin 202 


( 6 ) 
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The second line in Eq. ® gives the form of Jcp in terms of the parametrization Eq. (CD- 
Hence, if one of the entries hi 2 , ^13 or h 2 z vanishes, there is no CP violation in oscillation 
experiments. 

It is helpful to investigate the approximate form of h as implied by current data. We will see 
below that \ cp. 6 c — 0.22 is a useful parameter to describe the order of magnitude of the 
entries in niy. Since it holds to a good precision that 11^2 = ■\/l/2 (1 — A) and R = D 
(where D is of order one) [21], we can write in case of normal hierarchy ~ 
m\ ~ Am\D and in case of inverted hierarchy mf ~ Am\, rri^ ~ Am\ (1 + P A^). 
Ignoring the smallest mass state and the CP phase and setting for instance Ues = A 
and P ^3 = ■\/l/2 (1 — P A^), where A, B are real parameters of order one, we have in case 
of a normal hierarchy: 

/ A 2 A^ A 2 \ 

h ~ Ami ( ' ^ ^ I ’ 

and for the inverted hierarchy: 

/I A^ A^ \ 

h ~ Ami ( ' ^ ^ j ’ 

where the entries are to be understood as order-of-magnitude-wise. These two simple 
matrices are a very helpful guiding line in the search for promising candidates. 

Turning to the quasi-degenerate scheme, and in case of a normal ordering, we can write 
m 3 = uIq, m\ = mQ (1 — 77(1 — D A^)) and m\ = m^ (1 — 77 ), where 77 = Am\/rn^. Then, 
with the above dehnition of the mixing parameters, we hnd 

( 1 + 77 77 A^ 77 A^ \ 

1+77 77 + A^ j . (9) 

1 + 77 + A^ / 

A very similar expression is found for the inverted mass ordering. Note that the entries of 
order “1” in the above matrix are in our approximation exactly 1. This means, that the 
parameters in rriy are such that the leading terms in the diagonal entries of h are identical. 

We end our general considerations about Dirac mass matrices with the prospects of distin¬ 
guishing Dirac from Majorana neutrinos. Unfortunately, proving that neutrinos are Dirac 
particles is much more difficult than proving them to be Majorana particles. The only 
known process sensitive to the lepton number violation characteristic for light Majorana 
neutrinos is neutrinoless double beta decay [Sj. This process is triggered by the so-called 
effective mass (m), the ee entry of the Majorana mass matrix in the charged lepton mass 
basis (for details, see, e,g, the recent reviews [211221) • 

Suppose hrst that (except for the oscillation results) there is no signal indicating neutrino 
mass from direct laboratory searches or cosmological observations, i.e., we do not know 
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the exact neutrino mass scale. If the mass ordering is shown to be inverted, we know that 
— if neutrinos are Majorana particles — there is a lower limit on the effective mass 
of (m) > \//\m\ cos^ 6*13 cos 26 *i 2 — 0.01 eV. Hence, pushing the limit on the effective 
mass lower than this value would show that neutrinos are Dirac particles. In case of the 
normal ordering however, the effective mass can vanish even though neutrinos are Majorana 
particles. This would indicate that the ee entry of the neutrino (Majorana) mass matrix 
in the charged lepton mass basis is zero, which would be a very interesting observation 
by itself. Nevertheless, we have no handle on determining the nature of the neutrino. A 
possible way out of this dilemma would be to look for analogous processes of neutrinoless 
double beta decay. However, the associated branching ratios or cross sections are many 
orders of magnitude below current and future experimental bounds j2El- 
If we have however a dehnite value of the neutrino mass scale, be it from cosmology or 
direct laboratory searches, and do not hud a corresponding signal in neutrinoless double 
beta decay experiments, we could conclude that neutrinos are Dirac particles. To be more 
quantitative, if the neutrino mass scale is given by mo, the effective mass should be larger 
than roughly mo cos 26 *i 2 if neutrinos are Majorana particles. 

3 Minimal Dirac Mass Matrices 

In the remainder of the paper we shall consider minimal Dirac mass matrices, i.e. matrices 
possessing a certain number of zero entries. In the numerical search for successful matrices, 
we required the neutrino mixing observables sin^ 6*12,23,13 and R = Amg/Am^ to lie within 
their 2 a ranges given in Eq. O- The non-vanishing entries are varied within two orders 
of magnitude, i.e., we allow for up to a factor of 100 between them. We checked that the 
found textures are stable with respect to small perturbations, in the sense that the results 
stay unchanged as long as the “zero entries “ are one order of magnitude smaller than the 
non-zero entries. 

3.1 General Considerations 

It is fairly easy to see that the possible Dirac nature of neutrinos allows for at least as 
many zeros as in case of Majorana neutrinos. Consider in Eq. © the case V = \. Then 
we have rriy = U and therefore 




Del mi Ue2m2 Uesms \ 

11^2^2 11^37713 

Urimi Ur2m2 Ur3m3 / 


( 10 ) 


In the case of an extreme normal (inverted) hierarchy we can set mi (m3) = 0 and thus 
have three zeros in the hrst (third) column. Further setting 11^3 to zero, which is compatible 
with current data, adds another zero entry to for the normal hierarchical case. We will 
indeed encounter successful texture zero matrices with such a structure in the following 
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Sections. 


One can group Dirac neutrino mass matrices into classes which lead to the same matrix h. 
Matrices contained in these classes can be transformed into one another by a permutation 
of the right-handed neutrino helds, which corresponds to an exchange of columns in rriy^ 
i.e., 

rriy ml = m^Pi . ( 11 ) 

This leaves the rest of the Lagrangian invariant and leads to the same matrix h because 

h = m^ml h' = m^m^J = myPjP/mJ, = h . 

The orthogonal permutation matrices P* are: 

/100\ /010\ /00l\ 

Pi = 0 1 0 , P2 = 1 0 0 , P3 = 0 1 0 , 

\001/ \001/ \100/ 

( 12 ) 

/100\ /00l\ /010\ 

P4 = 0 0 1 , P5 = 1 0 0 , Pe = 0 0 1 . 

\010/ \010/ \100/ 

Note that we cannot put two matrices which are transformed into one another by a trans¬ 
formation of the left-handed lepton helds in one class, since we have already hxed the 
basis in which the mass matrix of the charged leptons is diagonal and therefore such a 
transformation would not leave the rest of the Lagrangian invariant. 

Finally, before proceeding with the search for successful texture zero mass matrices, it 
is worth commenting on the RG running of a Dirac neutrino mass matrix. The non¬ 
diagonal entries in the beta function for the Yukawa couplings are (in the charged lepton 
mass basis) proportional to [ 2 ZI where v is the weak scale. In principle, this 

term could destroy the texture zeros under consideration. However, the smallness of the 
Yukawa couplings, which are in the SM at most of order 10 “^^ and in the MSSM of 
order 10 “^^ (1 -|- tan^/ 3 ), renders the corrections negligible. Thus, ignoring possible SUSY 
threshold corrections, the stability of the texture zeros, which will be discussed in the next 
Section, is to a good precision guaranteed. 

3.2 Five Texture Zeros 

When the current neutrino data is taken into account, the maximal number of zeros in a 
Dirac mass matrix turns out to be 5 . There are in total 126 possibilities to put 5 zeros 
in m^, 18 of which allow for a successful explanation of the neutrino data. Their common 
feature is that their determinant vanishes and therefore they possess one vanishing mass 
eigenvalue (mi = 0 for NH and m3 = 0 for IH). Furthermore, there is no CP violation. 
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These 18 matrices fall into three classes which we denote A, B and B. Matrices belonging 
to class A predict the inverted hierarchy and t/gs = 0 whereas those belonging to B and B 
can accommodate both hierarchies and lead in general to non-vanishing Ue^. 

Class A is represented by the matrix: 


(13) 


which can be transformed into five other matrices with one of the permutation matrices Pi 
of Eq. (fT^ : 



APo = 


b a 0 

d 0 0 I , A P3 = 

e 0 0; 

6 0 a 
A Pi, = \ d 0 0 
; 0 0 


0 6a 

0 d 0 I , A P4 = 
0 e 0 


AP. = 




Next, we show that class A works only in case of an inverted hierarchy with the additional 
constraint of Ues = 0. Note that the matrices of class A predict for the resulting h a 
vanishing (23)-subdeterminant, independent of the choice of the variables a, 6, d, e. Using 
the dehnition h = U and requiring that the (23)-subdeterminant vanishes, 

leads to the condition (m^ (m^ 0^2 + '^ 12 ) + ^2 ~ 0- ^ normal hierarchy 

(i.e., mi = 0) this gives ?^ 2 ^ 3 Cg 2 Ci 3 = 0, which certainly cannot accommodate the data, 
since it implies Uei = 0. In the case of an inverted hierarchy (i.e., m 3 = 0) we end up 
with the equation m^m^s^g = 0 , which implies ^13 = 0 and leaves the other mixing angles 
unconstrained. 

In general, the formulae for the neutrino mixing observables in terms of a,..., Z are quite 
complicated, but in case of 5 zero entries one can obtain simple expressions. Taking a, 6, d, e 
real and positive we hnd 


tan^ 023 — ^ , sin^ 6^12 — | (1 + 


m + m — d^ — 

w 


p_ 

Am^ 


w 


|(a^ + 6^ + — tc) ’ 

where w = {a? + 6 ^ + d^ + — Aa^{d? + e^) 


(14) 


Typically, a, e and d have (in units of ■>/ Am\) absolute values of order one and 6 is of order 
A^. With these estimates, one obtains the approximate form of h from Eq. (jH)). Note that 
sin^di 2 can give the value sin^di 2 — 0.3 quite naturally. However, cancellations of order 
one parameters lead to the smallness of R which induces a certain hne-tuning. This will be 
the case for most of the solutions which predict the inverted hierarchy. There is typically 
a large (up to two orders of magnitude) gap between the entries of m^, and this implies 







that the inverted hierarchy will — in a bottom-up approach — always be associated with 
some tuning. 


The 12 other matrices able to reproduce the neutrino data can be divided into two very 
similar classes, which we denote B and B. They generate — independent of a,b,d,e — a 
vanishing hi 2 and hi^, respectively. Their representants are 

/0a0\ /0a0\ 

5=1 b 0 0 j and B = I b d 0 \ . (15) 

\ d e 0 J \ e 0 0 J 

Due to the vanishing of hi2 (^ 13 ) in class B (5) Jcp is also zero, i.e. there is no CP 
violation. From these conditions one can obtain an interesting correlation between the 
neutrino mixing observables. Setting the 12 element of h = U f/^ to zero gives^ 


\Ue3 


1 R sin 26*12 cot 6*23 

2 1 =p 5 sin^ 6*12 




-R sin 29 
2 


12 


cot 6*23 , 


( 16 ) 


where the — sign is for the normal and the + for the inverted mass ordering. The smallness 
of R implies that the condition hi 2 = 0 gives to a good precision the same correlation for 
both mass orderings. A similar relation can be deduced from hi^ = 0 for class B: 


\Ue3 


1 R sin 26*12 tan 623 

2 1 =p 5 sin^ 6*12 




-R sin 26* 
2 


12 


tan 023 • 


(17) 


The (close-to-)maximal value of 623 indicates that both, hi 2 = 0 and hi^ = 0, give ap¬ 
proximately the same result for Ue 3 - Moreover, the above relations Eqs. (IWH) are valid 
also for a non-vanishing smallest neutrino mass. Consequently, we can apply them also 
in later Sections of this work, when hi 2 or hi^ vanish for matrices with a non-vanishing 
determinant. We remark here that one can show that no other element of h except for hi 2 
and hi 3 can vanish without leading to conflicts with the neutrino data. We will comment 
on unsuccessful matrices in Section ITHl 

Inserting the best-£t values R = 0.036, sin^ 0 i 2 = 0.3 and sin^023 = 0.5, gives for both 
hi 2 = 0 and his = 0 as well as for both mass orderings that \Ue 3 \ — 0.017 which is about 
two orders of magnitude below the current 2(j upper bound for |t/e 3 p. Its value increases 
(decreases) with increasing (decreasing) sin^ 623 for a zero 13 entry in h (and vice versa for 
zero ^ 12 ). We plot in Fig.^the correlation between sin^ 023 and iDgsP for class B in case of 
the normal hierarchy. Fig. |21 displays the same for matrix 5, showing the slightly different 
behavior for these two observables as discussed above and implied by Eqs. (I16I17|) . The 
respective plots for the inverted hierarchy look identical, as expected from the approximate 
formulae discussed in this Section. As can be seen, |t/e 3 p stays between 10“'^ and 10“^; in 
particular, for the best-£t value of sin^ 623 = 0.5 we have 10“^ ^ jf/ggl ^ 5 ■ 10“^. 

^We have assumed C13 ^ 0. 
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Now we investigate the dependence of the observables on the variables a, 6, d, e with 
a,b,d,e > 0 in class B. Similar calculations for the class B reveal similar results and 
are therefore omitted. In case of normal hierarchy, one hnds: 


20 _ ya +0 + 


tan^ 023 = 






^ el 


(f + or {It + e^) 


UeZ = 


a(a^ — iP' + + w) 

(4 6^ + a^{a? — 6^ + d^ — + wY + d^{o? + 6^ + d^ — + tc)^) 

n |(a^ + 6^ + d^ + — tc) 


1/2 


2 y ^ I ty I I Wy j 

+ 6 ^ + d^ + + ic) 

where tc = ■\/(a^ + 6^ + d^ + — 4 (d^e 2 + 0^(6^ + e^)) 


(18) 

Typically, the normal hierarchy is achieved for parameter values (again in units of \/Am\) 
of order one for b, e and of order A for a and d. Inserting those values in mj, yields the 
approximate form of h from Eq. (ED, where however the 13 entry vanishes. With the order 
of magnitude of a, b, d, e we can obtain from Eq. m the following approximate expressions: 


tan^ 023 — K 

e 


sin^ dio — 


a^(&^ + e^) 


dV + e"(d" + d") 


f/e3- 


ad b 


r 


j— 


d^e^ + dYb^ + eY ,2 

— wrrf —^ ’ 


(19) 


where we have dehned new order one parameters d = a/X and d = d/X. As can be seen, Uez 
and R are small because they are proportional to and, again, sin^ di 2 can be naturally 
of order one. 

In case of inverted hierarchy, we get: 


tan^ 02.3 = To , = 


de 


sin^ dio = i 


b^ ’ 

^4 ^ 2/^2 


R = 


y^a^ 6 ^”+”e 2 (a 2 ”+”d^ 
e^(a^ — 26^ + d?) — b‘^{a? — b^ — d?] 
( 6 ^ + e^w 

w 


w) ’ 


i(a^ + 6 ^ + d^ + 
where w = (a^ + 6^ + d^ + e^Y — A^dXe^ + a^ih^ + e^)) 


( 20 ) 


In order to explain the neutrino data one needs now that a, 6, e are of order one and d 
of order A^ (leading approximately to Eq. (jHD with a vanishing 13 entry). Analyzing Eq. 
CT shows that R owes its smallness to a cancellation of the order one coefficients which 
is similar to the matrices in class A discussed above. 

To conclude this Section, both the normal and the inverted hierarchy can be accommodated 
with Dirac neutrino mass matrices containing 5 zero entries. No observable CP violation 
is predicted and is either exactly zero (then the inverted hierarchy is present) or very 
small, corresponding to |t/e 3 p — 10“^ ... 10“^. A non-vanishing value of |f/e 3 p above 10“^ 
will therefore rule out any Dirac mass matrix with 5 zero entries, while a value between 
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10“^ and 10“^ excludes class A. Class A is also ruled out when the neutrino mass spectrum 
is found to be normally ordered. The three classes of matrices A, B and B are furthermore 
excluded when either leptonic CP violation is present or when a positive signal of neutrino 
mass in laboratory searches or by cosmological observations is found. Distinguishing B 
from B will be very difficult, unless 6^23 deviates sizably from vr/d. We give in Tabled the 
classes of matrices which are successful in explaining the neutrino data. Together with the 
number of members in the class we give their phenomenological consequences and an order 
of magnitude estimate of the non-zero entries. 

3.3 Four Texture Zeros 

As shown in the previous Subsection, Dirac mass matrices with hve zero entries can ac¬ 
commodate neutrino data, but predict in general vanishing CP violation. To take this 
effect into account enforces one to allow for an additional non-zero entry in m^. As for 
the case of 5 zero entries, there are 126 possible matrices, 69 of which are successful in 
explaining the data. Those 69 belong to 12 classes. Among those are 4 classes with a 
vanishing determinant. From those 4 there are 3 classes which can accommodate leptonic 
CP violation. Let us start with the CP conserving matrices. 

The hrst class, denoted C and containing 3 members, is represented by 

/ a b d\ 

C= e 0 0 . ( 21 ) 

\f 0 0 / 

The resulting structure of h is analogous to the one of class A in the previous Subsection, 
hence one massless neutrino, zero tAs and the inverted hierarchy are predicted. At least 
one of the entries has to be suppressed by in order to reproduce the approximate texture 
of Eq. (jHl) . 

Other classes have a vanishing 12 or 13 entry of h and therefore, as argued above, a 
vanishing Jarlskog invariant. The resulting correlations between the observables, as given 
in Eqs. fjl6ll7j) and shown in Figs. [Hand El hold again in these cases. Consider the following 
18 matrices, which fall into the three classes and predict hi 2 = 0 : 

/a 00 \ /a 00 \ / a b 0 \ 

Di= 060 ^ D2 = \ 0 b d \ ,D 3 = 00 d . ( 22 ) 

\d e f j Ve 0 f ) \e 0 f ) 

From Eq. (|TH|l . one has \Ue^\ ~ \R sin26^12 cot 6 ^ 23 , whose behavior is displayed in Fig. El 
To accommodate a normal hierarchical spectrum, the second class D 2 has typically a, b, e 
of order A and d, f of order one. If e is of order A^, d of order A and the other parameters 
of order 1, it turns out that the spectrum can become quasi-degenerate. This can be seen 
from the fact that the determinant of the matrix, det D 2 D\ = a? 6 ^ p is then close to the 
third power of the trace, Tr D 2 D\ = [a? + b"^ + d?‘ -\- -\- /^). Recall that in a bottom-up 
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approach, the case of quasi-degenerate neutrinos is hne-tuned, since the diagonal elements 
of h have to be identical. The inverted mass ordering is also possible, but as in the classes 
B and B discussed for matrices with 5 zeros, this is connected with a large spread of 
parameters. To be precise, the entries b and e of D 2 have to be of order and a, d, f of 
order 1. 

The following three classes h)i 2,3, giving in total 18 matrices, generate a zero in hi^: 

/ a 0 0 \ /a00\ / a b 0 \ 

Di = i b d e \ , D 2 i b d 0 \ ,i)3=|(i0ej . (23) 

Vo / 0 / Vo e // VO 0 f J 

Thus, for them Eq. (HID holds and the behavior of Ues and sin^ 623 can be seen in Fig. ^ 
Similar statements about the neutrino mass scheme as for the classes £* 1 , 2,3 can be made. 


We furthermore hnd two classes denoted E and E, which contain six matrices each. Rep- 
resentants look like: 

/ a b 0 \ / a 0 b \ 

E = I d 0 0 j and E = I d e 0 \ . (24) 

Ve 0 f J \f 0 Oj 


Their common feature is a non-vanishing determinant and Jcp 
one hnds that 

hi2 _ 

hl3 h23 

for class E while for class E it holds 


0. From the matrices 
(25) 


hl2 _ ^32 
hl3 h33 

Inserting Eq. (|2SD in the dehnition of Jcp from Eq. © gives for class E that 


(26) 


Jcp cc Im(/!i2ft23/i3i) ^ Im(/i 22 /!i3 ^ 31 ) ^ 0 , 


(27) 


since h is a hermitian matrix. A similar equation can be written down for class E. This 
explains why there is no CP violation for the classes E and E. 

The conditions resulting from Eqs. ( 1 ^ and (EED can be used to obtain an interesting 
correlation between the neutrino mixing observables. Starting with class E, the relation 
Eq. (I25D leads for the normal mass ordering to: 


|t^e3| = 


AmQ{Am\ + ml) sin 26 *i 2 tan023 


Am^Am^ + (2 Am\ — AmQ)mf + AmQ{Am\ + mj) cos20i2 


tan 023 tan 012 


for m? = 0 


(28) 


jR sin 2012 tan 023 for 
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where we gave the limits for the normal hierarchical and the quasi-degenerate spectrum. 
Note that the limit for m\ ^ A/tIq a, which corresponds to a quasi-degenerate spectrum, 
is just the correlation of the parameters Eq. 03 for a vanishing 13 element of h. The fact 
that the value of f/e3 in case of zero mi is too large indicates that the matrices of class E 
lead to a lower limit on the neutrino mass if the spectrum is normally ordered. Indeed, 
one can solve Eq. ( 0 H|) for mf/ArriQ, which gives: 

ml _ sin 26*12 sin 6*23 — 2 cos^ 6*12 cos 6*23 lUg^l 
AtUq 2{1 — R sin^6*i2) cos6*23 \ Ue^\ — R sin6*23 sin 26*12 

I for Uez = 0 

I cos6*12 ^sin 0 i|_^£^ _ for \Uez\ '> R 


Consequently, m\/Am?Q is minimized for the largest possible \Uez\, and is then given by 
ml ^ i Am^ sin 26 *i 2 tan6*23/]t/e3|- Putting in numbers gives mi ^ 0.01 eV. The mass 
spectrum can become degenerate when the denominator of Eq. (Ei goes to zero, which 
occurs when Uez takes the value given in Eq. (ED- This value is about 0 . 02 . Lower values 
of \Uf.z\ render m^/AmQ negative and are not possible. Hence, \Uf.z\ — \R sin26*12 tan6*23 
is the lowest possible value for class E and the normal mass ordering. 

In case of an inverted mass ordering Eq. ( 123 ) leads to 


|f^e3 


iAm^mg sin 26*12 tan 6*23 
Am^(m3 + Am^) + Aml^{Am\ — rril sin^ 6*12) 


r'^ 


0 


1 

2 


R sin 26*12 


tan 6*23 


for m| = 0 
for m\ ^ AmQ ^ 


(30) 


In this case there is no lower limit on the smallest mass state because \Uez\ goes to zero 
when m3 goes to zero. One can again solve Eq. dH for the ratio of the smallest neutrino 
mass m3 divided by the solar mass squared difference, which gives 


ml ^ _ 2(1 - R) cos 6*23 1^631 _ 

Am^ R (2(1 + R sin^ 6*12) cos6*23 \Uez\ — R sin 6*23 sin 26*12) 




0 for Uez = 0 
for > R 


(31) 


The denominator vanishes when Uez fulhlls Eq. ca, thereby corresponding to quasi¬ 
degenerate neutrinos. The expression for ml/Am^ becomes negative when Uez is larger 
than roughly 0 . 02 . Hence, | 17 e 3 | — ^R sin 26*12 tan 6*23 is the largest possible value for class 
E and the inverted mass ordering. 
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We plot in Fig.ISlfor both mass orderings and class E a scatter plot of \Uez\ against the ratio 
of the smallest neutrino mass squared divided by AtUq. The lower limit of | 17 e 3 |, roughly 
0.02 for the normal mass ordering, occurs for large mf/Amg, i.e., for quasi-degenerate neu¬ 
trinos. It is given, as discussed above, by R /2 sin 26^12 tan 6*23. The qualitative discussions 
given above are nicely conhrmed by the plot. 

For the matrices belonging to class E we have the relation hi2/hi^ — h^2/hzz = 0 , which 
leads to conditions for the observables very similar to Eqs. They are obtained 

from those equations by replacing tan 6^23 with cot 6^23. 

The remaining classes of 4 zero mass matrices predict non-zero Ue3 and non-vanishing 
CP violation, which — from the phenomenological point of view — are perhaps the most 
interesting ones. There are three classes, all of which have a vanishing determinant, i.e., 
one neutrino is massless: 

/a00\ / a b 0 \ / a b 0 \ 

Fi=|6(ioj ,F2=|dOoj ,F3=|(ieoj . ( 32 ) 

Ve / 0/ \e f 0/ \f 0 0/ 

Unfortunately, there is in general no interesting correlation between the parameters and 
all observables, including the CP phase, can take any value inside their allowed regimes. 
Again, the reader can easily obtain the rough order of magnitude of the entries in Ti^2,3 in 
order to reproduce the approximate textures given in Eqs. (PHD. However, in contrast to 
the case of 5 zero entries, the position of the large (i.e., order 1) and small (i.e., order A or 
A^) terms is not unique and there is some freedom. Examples are given in Tables |21 and |Hl 
together with the implied phenomenology and the number of matrices in a given class. 

We encounter here the simple example given in the beginning of Section 13.11 If the pa¬ 
rameters in class Ei are such that eI Ei is diagonal, we have V = 1 and can accommodate 
the normal hierarchy together with Uez = 0 , see the discussion after Eq. m- This would 
require in Eq. that db* + f e* = 0 . 

We can rule out class C with a non-vanishing value of 11^3 or when the mass ordering 
is normal or when laboratory searches or cosmological observations hnd a signal corre¬ 
sponding to a non-zero neutrino mass. Classes £*1,2,3 and £*1,2,3 (which are difficult to 
distinguish) are ruled out for a value of Ues larger or smaller than ~ P /2 sin 26*12. They 
can in principle be distinguished from the 5 zero classes B and B because they allow for 
a non-vanishing smallest neutrino mass. Classes E and E are difficult to distinguish and 
are excluded when the mass ordering is normal (inverted) and Ues is smaller (larger) than 
P /2 sin 26 * 12 . All the classes C, £*1,2,3, -^1,2,3, E and E are ruled out if there is leptonic CP 
violation. Finally, classes £1,2,3 are ruled out if all three neutrinos are found to be massive. 
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3.4 Three and less Texture Zeros 


As we have seen, already for matrices with 4 zeros the results for the observables can be 
somewhat arbitrary (see the discussion for the classes -^ 1 , 2,3 after Eq. dSl). The major 
part of the successful matrices with 3 or less zeros falls into this category. Counting the 
number of possibilities for three zero entries, one arrives at 84 of which only 6 generate one 
vanishing eigenvalue. Three of these can accommodate all data and are represented by: 

/ a b 0 \ 

d e 0 . (33) 

\f 9 Oj 

No correlation of observables, except that mi or m 3 vanishes, is predicted. This provides 
a possibility to exclude them. Again, if the parameters in this matrix are such that mj, 
is diagonal, we have V = 1 and can accommodate both hierarchies, compare with the 
discussion after Eq. m- 

The remaining 78 matrices without vanishing eigenvalues can be divided into sixteen 
classes. Two of those, each containing three matrices, fulhll the relation hi 2 = 0 and 

/ 00 a\ / 0 a b\ 

I 6 d 0 j and j d 0 0 j . 

Ve f 9 J \e f g) 

They lead to the same results as the matrices belonging to classes 2,3 (hence are not 
distinguishable from the latter), in particular they predict CP conservation, and have the 
same correlation between Ues and sin 623 as given in Eq. (ITT)|) and Fig. |21 There also exist 
two classes with his = 0 whose representants are: 

/ 00 a\ / 0 a b \ 

I 6 d e and | d e / . 

\f 9 0 J \g 0 0 J 


They lead to the same results as the matrices belonging to classes £* 1 , 2,3 and consequently 
are not distinguishable from them and only very hard to distinguish from the above two 
classes. The presence of CP violation will exclude the last four classes. All the other 
classes contain six matrices each and are represented by: 


Gi 

C5 


0 0 a \ 

0 6 d j 

e f 9 J 

0 a b \ 
d 0 e j 

0 / gj 


G2 = 

Cq = 

Cg = 


00a\ / 0 a b \ 

b d ej,G3=|o 0 dj,G4 

^ f 9 / \e f 9 / 

0 a b \ / 0 a b \ 

d 0 e , Gr = d e / , Gs 

f 0 g J \0 0 g J 

a b d \ / 0 a b \ 

0 6 / 1 , Gio = I d 0 e j . 

0 0 g J \f 9 0 J 


0 a b \ 

0 d e j , 

/ 0 g J 

a b d \ 

0 0 e , 

0 / gJ 


(34) 
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None of these reveal simple correlations among the observables. 

Similarly, the 36 possible matrices with two texture zeros and the nine ones with only 
one zero (of course, all of them are allowed) do not show any relations among the mixing 
angles, the CP-phase and the ratio of the mass squared differences. 

3.5 Symmetric Matrices with Texture Zeros 

The discussion so far assumed a general structure for the Dirac neutrino mass matrices. In a 
given model based on some symmetry, however, symmetric matrices might arise. Extending 
our discussion to the case of symmetric Dirac mass matrices is rather straightforward. It is 
known from the discussion of Ref. na that Majorana mass matrices, which are necessarily 
symmetric, can only have two independent zero entries. It turns out now that if the three 
zero Dirac mass matrices in the classes G-i (see Eq. (ED) are symmetric, they can reproduce 
six of the seven allowed two zero Majorana mass matrices from Refs. na. In the language 
of those references, the two zero Majorana mass matrix Ai corresponds to the symmetric 
form of matrix Gi from Eq. (IH^ . Analogously, it is easy to see that 6 * 2 , Gs, G-j from Eq. 
dMD correspond to A 2 ,B 3 , from Refs. na, respectively. The matrices belonging to class 
Gio reproduce cases Pi and P 2 . We display this correspondence in Table El Case G from 
na would correspond to a symmetric mass matrix with two zero entries, whose discussion 
we omitted in this work, since in general no interesting correlation of variables arises from 
this possibility. The phenomenological consequences of two zero Majorana mass matrices 
are rather well known na and we have nothing to add here and refer to those works for 
details of their phenomenological implications. Nevertheless, we included the approximate 
formulae out of which the phenomenology follows in the Table. 

3.6 Unsuccessful Matrices 

We hnish this Section with a few comments on matrices which are unsuccessful in explaining 
the neutrino data. 

In case of 5 zeros in rrty, there are three categories of matrices which do not work. First, 
there are 90 matrices which lead — independently of the choice of the four complex pa¬ 
rameters — to the vanishing of two mixing angles in U. In the second category, there are 
twelve matrices which generate a zero ( 12 )- or (13)-subdeterminant in h. A vanishing of 
the ( 12 )-subdeterminant leads for the normal hierarchy (mi = 0 ) to Uri = 0 and for the 
inverted hierarchy (m 3 = 0 ) to Urs = 0 . Analogously, if the (13)-subdeterminant vanishes, 
one ends up with = 0 for the normal mass ordering and = 0 for the inverted 
ordering. All these conditions are of course not compatible with the data. In the last 
category, there are six matrices predicting ^23 = 0. Phenomenologically ^23 = 0 leads to a 
too large value for Pe 3 for the normal as well as the inverted hierarchy. Hence, 108 out of 
126 possible matrices with 5 zeros are unsuccessful in explaining the neutrino data. 

For matrices with four zeros we hnd that about half of the 126 possible matrices, i.e., 57 
matrices, do not work. 27 matrices — 9 of which have a non-vanishing determinant - 
are not successful since two of the mixing angles are zero. Furthermore, three matrices 
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generate a vanishing (12)-subdeterniinant in h and further three ones a vanishing (13)- 
subdeterminant. Apart from these there are 18 matrices leading to h 23 = 0. Finally, there 
are six matrices fulhlling the equation hii/hi 2 = h^i/h^ 2 - In this case \Ue^\ turns out to be 
always bounded from below by ~ 0.3 and hence one cannot accommodate the data. 

In Table El we give typical examples of matrices with 4 and 5 zero entries which are not 
able to explain the data. 

Only 9 of the possible 84 three zero textures do not work. Three of them enforce two 
mixing angles to be zero and the six other ones generate ^23 = 0 . 


4 Realization of the Textures 

In this Section we indicate models which might give rise to the textures found to be able 
to reproduce the data. They originate either from the Froggatt-Nielsen mechanism j2Hj — 
explaining the hierarchy among the non-zero matrix elements — or from more complicated 
flavor symmetries. As an example, we choose the global non-Abelian discrete symmetry 
Di X Z 2 which is broken at the electroweak scale by enlarging the Higgs sector through 
scalars transforming non-trivially under the horizontal symmetry^. A global continuous 
symmetry broken at high energies is assumed to be the origin of this discrete symmetry. 
Rather than working out the details of the models, we merely reproduce the structure of 
the mass matrices, so that this Section should be understood as a proof of principle for the 
Dirac neutrino mass textures under consideration. One still faces the problem of explaining 
the relative smallness of neutrino masses with respect to the charged lepton masses. To 
achieve this, we have to rely on some tuning of the parameters of the models based on 
the discrete symmetry. In the model based on the Froggatt-Nielsen mechanism, we can 
explain the relative smallness by appropriate choice of charges and by working with two 
scalar fields, whose vacuum expectation values display a hierarchy. 

4.1 Mass Textures from the Non—Abelian Discrete Symmetry D 4 

X Z 2 

The group D 4 has already been discussed as a possible flavor symmetry by several authors 
[SDunuisainni. Other dihedral and two-valued dihedral groups have also been considered 
in the literature jHl] . Here we discuss D 4 for the first time in the context of Dirac neutrinos. 
Mathematical details of D 4 are delegated to the Appendix. 

4.1.1 Obtaining Class A from D 4 x Z 2 

Consider the Dirac neutrino mass matrix to be of the form of Eq. m and the charged 
lepton mass matrix to be diagonal. As shown in the following, this setup can be realized in 
the framework of the Standard Model extended by the non-Abelian flavor symmetry D 4 

^Of course there are several ways to realize texture zeros in mass matrices, see, for instance |29| . 
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X Z^- The left-handed leptons, right-handed charged leptons and nentrinos are assigned 
to be: 


1 + o-t 

-V- 




tr 


ll,2- 


{^ v ) r -, {^ t ) r } ~ 




The left-handed lepton donblets are Le = ^ ^ and analogonsly for and The 

Higgs sector of the model contains 5 SM-like SU{2)l donblets, two of which form a two- 
dimensional representation nnder ^4 and the three other ones a one-dimensional one: 





where ij denotes the representation of dimension i of D 4 and the snperscript ^ the trans¬ 
formation property nnder Z 2 . Therefore, a held (j)'^ transforms nnder Z 2 as whereas 
<p~ transforms as — 

Using real representation matrices for the two-dimensional representation, this leads to 
the following tree level mass matrices 

/ (0^) K2 (0^) 0 \ 

nil, = I 0 K 3 (' 0 ^) 0 and 

\ 0 0 / 

/ K4(0r) 0 0 

mi= \ 0 K 5 (0r) +(0r) 0 

\ 0 0 Kg (0]") - Kg (0r) 



which corresponds to the Dirac mass matrix of class A with a diagonal charged lepton mass 
matrix. The Ynkawa conplings are denoted with Ki and the Higgs helds 0^ and ipf have 
acqnired their vacnnm expectation valnes (vevs). Obtaining a hierarchy for the entries 
of the mass matrices reqnires for this and the next two examples, in general a detailed 
analysis of the Higgs potential, which is beyond the scope of the present work. Still one 
can give some order-of-magnitnde estimates for the vevs and Ynkawa conplings. This 
illnstrates the tnning reqnired to explain the relative smallness of nentrino and charged 
lepton masses. Invoking the framework in a scenario involving higher dimensions might 
explain the overall snppression of nentrino masses along the lines of u u This is well 
beyond the scope of onr work, which is snpposed to deal mainly with the phenomenological 
conseqnences of textnre zeros. Let ns give an example of the reqnired valnes of conplings 
and vevs: with v being the electroweak scale one can choose ( 0 ^) = 0 ( 10 “^^) v , ( 0 ^) = 
0 ( 10 “^^) n,( 0 ^) = 0 ( 10 “^^) V, ( 0 ]“) = 0 ( 10 “^) n, ( 0 J) = 0 ( 10 “^) v for the vevs and 
Ki = 0(10“^°), K 4 = 0(10“^), ^ 2 , K 3 , Kg = C>(1) and Kg = 0{1). Note that with these 
choices there has to exist one fnrther Higgs donblet whose vev is of the order of v. For 
example, one might introdnce a Higgs held transforming as which cannot conple to the 
leptons becanse of their assignments nnder x Z 2 , bnt can give adeqnate masses to the 
qnarks. 
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4.1.2 Obtaining Class B from x Z 2 


To obtain class B from the group ZI 4 x Z 2 we choose the following assignments for the 
leptons: 



while for the Higgs helds we now require 6 doublets: 



These choices lead to the following mass matrices: 

/ 0 K2 {(pt) - Ks {(pt) 0 

= ^2 ( 02 ) + ^3 (03 ) 0 0 1 and 

\ «i( 0 ri «i( 0 ^) ( 

K 4 (0r) + ^^5 (04 ) 0 0 

mi= \ 0 «4 (0r) - (0r) 0 

0 0 Kq ( 01 ") 


The hierarchy among the elements in the mass matrices can be reached in a similar way 
as explained in the previous example. Note that here we have different vevs appearing in 
the matrices governing rriiy and mg, which slightly simplihes the issue. 


4.1.3 Obtaining Class Di from ZI 4 x Z 2 


A texture belonging to class Dg can be reproduced by the following assignments for the 
leptons which is very similar to the assignment used to create textures belonging to the 


class B\ 



For the Higgs particles we choose: 



yielding the mass matrices 


m,, = 


mg 


1^1 {(Pt) + 1^2 {(pt) 0 \ 

0 Hi {(pt) - 1^2 {(pi) 0 j and 

Ksilpt) Hsi^Pt) «4(0r) / 

/ Kg, {(pp) + Kq (0J) 0 \ 

0 Kg {(pp) - Kq {(p-^) 0 

\ 0 0 Kt{(PP) J 
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Again we can reproduce the hierarchy among the elements and the relative smallness of 
the neutrino masses with an interplay of vevs and Yukawa couplings. 


4.2 Mass Textures from the Proggatt—Nielsen mechanism 

In this Subsection we show how one can generate the textures under consideration in 
a supersymmetric framework by assuming the existence of two global horizontal flavor 
symmetries denoted U{1)a and U{1)b- As an example, we generate class A of Section 
EH In addition to the two global U{1) symmetries, we require the presence of the cyclic 
symmetry Z 3 acting only on the charged leptons to preserve their mass matrix to be 
diagonal. Under this Z 3 symmetry, the right- and left-handed charged leptons transform 
according to ~ where 00 = . The U{1)a,b are both broken at high-energy 

scales Ma and Mb, respectively. There are two scalar helds a and p, which are gauge group 
singlets. The held a (p) carries charge —1 under U{1)a {U{1)b) and, by acquiring a vev, 
whose magnitude is of the order Ma {Mb), breaks the horizontal symmetry spontaneously. 
Small parameters associated to this framework are given by the ratio of those vevs with 
some high scale M, which corresponds to the mass of the additional heavy vector-like 
fermions, which are integrated out. The t/(l)-charges are chosen as follows: 

{Le, L^, Lt} ~ {(3^, Os), (5a, Os), (5a, Os)} , 

{^R, tr} ~ {(8a, Os), (7a, Os), (5a, Os)} , 

{{^e)R, {j^tJ,)R, {^t)r} ~ {(3a, 6 s), (5a, 6 s), (Oa, 6 s)} . 

The MSSM Higgs doublets and are assumed to have zero charges under both U{1) 
symmetries. For the helds p and a we choose the following possibilities for their charge 
assignment and their vevs: 


a ~ (—1a, Os) with {o')/M = A ~ 0.22 , 

P ~ (Oa, — Ifi) with {p)/M = e ~ 10“^ . 

The relative smallness of {p)/M is required to explain the hierarchy between the charged 
lepton and neutrino masses and between the leptons and quarks, whereas {cr)/M is respon¬ 
sible for the mass hierarchy between the generations of the charged leptons and neutrinos, 
respectively. We end up with: 


m. 


= 0{e^ 


/ 0(1) 

0(A2) 

/ 

( 0(A5) 

0 

0 \ 

0 

0(1) 

0 {Hu) and me = 0(1) 

0 

0(A2) 

0 

V 0 

0(1) 

0 \ 

V 0 

0 

0(1) / 


(H,) 


Note that the zeros in my are preserved from being hlled, since we are working in a SUSY 
framework, in which the conjugated held of a, which would carry the charges (1 a, Os), is 
absent (“SUSY zeros”). Finally we would like to note that the relative smallness of the 
neutrino masses has its origin in an appropriate choice of the charges and the small value 
of {p)/M with respect to {a)/M. 
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5 Conclusions, Summary and final Remarks 

Apart from theoretical prejudices, we have no incontrovertible evidence for the Majorana 
nature of neutrinos. If they are Dirac particles, just like all other known fermions, one 
might be interested in the form of the Dirac mass matrix implied by current data. In 
particular, the minimal allowed structure of might be of special interest. Here, the 
situation is analyzed in terms of putting as many zeros in rrii, as possible. As a proof 
of principle, we have also outlined possibilities to obtain the successful Dirac neutrino 
mass matrix textures from models based on either discrete non-Abelian symmetries or 
on a Froggatt-Nielsen scenario. We found that up to 5 zero entries in m,^ are allowed 
by current data. All of those candidates imply CP conservation, one massless neutrino 
and their phenomenology is summarized in Table d Most successful matrices with 4 zero 
entries also predict CP conservation, their physical implications are summarized in Tables 
El and 01 Apart from a few exceptions, which then reproduce the phenomenology of already 
discussed matrices with 5 or 4 zeros, all matrices having 3 or less zeros imply no correlation 
between the neutrino mixing observables. Typically the inverted hierarchy requires that 
there is a spread between the non-zero entries in ruj, of two orders of magnitude, whereas 
the normal hierarchy can be successfully obtained with a rather mild hierarchy. We would 
like to remark that some matrices predict ^13 = 0 , but none of them implies in general 
^23 = 7 r/ 4 . This would require typically that two of the non-zero entries in rUi, are identical. 
It is hnally interesting to note that there are basically only three signihcantly different 
predictions coming from Dirac mass matrices with zero entries. Those are the inverted 
hierarchy with zero fAs and the two possibilities given in Eqs. (nmTi) and (12813011 . Note 
that once there is an interesting correlation between observables, CP is automatically 
conserved. The characteristic value of Ue 3 which can play the role of a discriminator 
between those possibilities, is given by iLAsI ~ |i?sin 26 * 12 . 
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Class 

N 

Condition 

Phenomenology 

Order of 
NH 

magnitude 

IH 


f a b 0 ^ 
0 d 0 
e 0 ) 


6 

^22 ^23 _ Q 

^32 ^33 

only IH 

Ue3 = 0, JcP = 0 

— 

1 

(1 A^ 0^, 
0 1 0 
Vo 1 0 ) 


B=\ 

/ 0 a 0\ 

boo 
d e 0 J 


6 

hi2 = 0 

Jcp = 0 

Ue 3 ^ f sin 26 'i 2 cot 623 
Eq. (HH), Fig. El 

1 

(0 A O '] 
1 0 0 
Vl A 0) 


1 

(0 1 ov 
1 0 0 
Vl A 2 0 ) 



/ 0 a 0^ 

b d 0 

iy e 0 0 y 


6 

hl3 = 0 

Jcp = 0 

Ue 3 ^ f sin 29i2 tan 623 

Eq. dni), Fig.n 

1 

fo A 0] 
1 A 0 

Vi 0 0 ) 


1 

(0 1 oV 

1 A^ 0 

Vi 0 0 ) 



Table 1: Successful matrices with five zero entries. Shown are a typical representants of the 
class, the number N of distinct matrices in that class, the characteristic condition giving 
its phenomenology and an order of magnitude estimate of its entries. We also indicate from 
which equation the phenomenology results. NH denotes normal hierarchy and IH inverted 
hierarchy. Note that the smallest mass state is always zero because the determinant of m„ 
vanishes. 


Class 

N 

Condition 

Phenomenology 

Order of 
NH 

magnitude 

IH 


^ a b d \ 
e 0 0 

Kf 0 0/ 


3 

^22 ^23 _ Q 

^32 h33 

identical to A 

— 

( 

1 A2^ 

1 0 0 
VI 0 Oy 


F .= 

/ a 0 0 V 
b d 0 

U / oy 


6 

— 

Jcp 7^ 0 

f/e37^0 


/A 0 0\ 



(1 0 0\ 
A2 1 0 
Va^ 1 0 ) 


F2 = 

f a b 0 ] 
d 0 0 

f oy 


6 

— 

same as above 

1 

A 0^ 
1 0 0 
Vl A Oy 


( 

(\^ 1 0^ 
1 0 0 
Vl A^ Oy 


F3 = 

f a b 0 ] 
d e 0 

\f 0 OJ 


6 

— 

same as above 

1 

A 0^ 
1 A 0 
Vl 0 Oy 


( 

(\^ 1 0^ 
1 A^ 0 
VI 0 Oy 



Table 2: Same as Table ^ for matrices with four zero entries and one vanishing mass 
eigenvalue. 
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Class 


N 


Condition 


Phenomenology 


Order of magnitude 


NH 


IH 





hio — 0 


similar to B 

Ue3 ~ f sin 2012 cot 023 

Eq. (Uni), Fig. El 



Do = 



hi2 — 0 


same as above 



^3 = 



0-12 — 0 


same as above 




D,= 



hi3 — 0 


similar to B 

Ue3 - f sin 2012 tan 623 

Eq. dni), Fig.^ 



Do = 



hi3 — 0 


same as above 



D 3 = 



h^3 — 0 


same as above 



NH: mi ^ 0.01 eV and 
Ue3 > f sin 2012 tan 623 

Eqs. Fig.El 


E = 



0-12 

0-13 


Ihi 

023 



IH; masses free and 
Ue3 < f sin 2012 tan 623 
Eqs. 


, Fig. H 


E = 



0-12 

0^ 


032 

0^ 


as class E with 
023 —> 023 + 7r/4 


(X^ 

0 

A\ 


0 

1\ 

1 

A2 

0 

1 

A2 

0 

v 

0 

0/ 

v 

0 

0/ 


Table 3: Same as Table El for matrices with four zero entries and no vanishing mass eigen¬ 
value. All classes predict Jcp = 0. 
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Symmetric Class from Eq. (14411 

Results of [T 7 ] 

Order of magnitude 


^ 0 0 a \ 
0 b d 
y a d e J 


Case Ai 
only NH 

ITT 1 VR Q sin 26 »i 2 

|Ce 3 | - 2 cot (^ 23 ^/^;^^^^ 


/O 0 A\ 

i-U 


G2= I 

^ 0 a 0 \ 
a b d 
y 0 d e J 


Case y42 
only NH 

lUesl^ ^tan 023 

1 2 Vcos 201 <;> 


/O A 0 \ 

(0;;) 


G 3 = I 

f a 0 b \ 
0 0 d 
y b d e J 


Case B3 

QD; both orderings 
iDe 3 i 

1 1 cosd ll+tan^0io 

1 

(\ 0 A2^ 

0 0 1 
VA^ 1 Ay 


Gj = j 

/ a b 0 \ 
b d e 
y 0 e 0 / 


Case B4 

QD; both orderings 

ITT 1 PI tan 2023 1 tan0i2 

1 e 3 | 1 cos (5 Il+tan^0i9 

1 

(\ A^ 0^ 
A^ A 1 
vO 1 Oy 


Gio = 

and 

/ a 0 b ^ 
Ode 
^ b e 0 J 
/ a b 0 ^ 
\ b 0 d 
\ 0 d e ^ 


Case B2 

same as B4 (or G7) 


^ 1 0 A^^ 

0 A 1 
^A^ 1 0; 


Case Bi 

same as R3 (or G3) 


H A^^ 

A^ 0 1 

^ 0 1 Ay 


/ a b d \ 

1 6 0 e 1 
\ d e 0 J 

Case C 

QD; both orderings 

\TT 1 cot 2012 cot 2023 

___cos _ 

( 

^ A^ A' 
A^ 0 1 



Table 4 : Symmetric Dirac mass matrices from Eq. (ED together with the corresponding 
two zero Majorana mass matrices from Refs. [TZj. Class C from El would correspond to 
a symmetric Dirac mass matrix with two zero entries. 
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Class 

N 

Condition 

Phenomenology 

matrices with one vanishing eigenvalue 


/ 0 6 0 \ 
0 d 0 
\ a e 0 j 


6 


hii hi2 
^21 ^22 

= 0 

NH; Uri = 0 

IH: Ur3 = 0 


f a b 0 \ 

d 0 0 

0 e 0 yi 


6 

^23 — 0 

Ces too large 

1 

^ a 0 0 ^ 
b d e 

{f 0 OJ 


3 


hii hi 3 
hsi ^33 

= 0 

NH; U.i = 0 

IH: = 0 

matrices with all 7^ 0 

1 

( a b d \ 
0 0 e 

VO / oy 


6 

^23 — 0 

similar to second matrix 
in this Table 

1 

a 0 0 V 
b Q d 

v^ / oy 


6 

hii _ hsi 
hl2 ^32 

t/e3 too large 


Table 5 : Typical matrices with four and five zeros unable to accommodate the data. For 
discussion see Section ESI 



























A Group theory of i^4 


The group theory is given by the following character table, generator relations and matrices 

ISSj: 



classes 


Cl 

C2 

C3 

C4 

C5 

G 

1 

A^ 

B 

AB 

A 

°C- 

1 

1 

2 

2 

2 

°hc. 

1 

2 

2 

2 

4 

ll 

1 

1 

1 

1 

1 

I2 

1 

1 

-1 

1 

-1 

I3 

1 

1 

-1 

-1 

1 

I4 

1 

1 

1 

-1 

-1 

2 

2 

-2 

0 

0 

0 


Here we denoted with Ip 2 3 4 2 the irreducible representations of G are the 

generators of the group, °di is the order of the class C, and °h a is the order of the 
elements in class Cj, where i = 1 ,..., 5 . The generators A and B for 2 are 

A=(J and . 

and fulhll the generator relations 

= 1 , = 1 and ABA = B . 


From the character table one can read off the Kronecker products: 


X 

Ip I2 I3 I4 

ll 

Ip I2 I3 I4 

I2 

I2 ll I4 I3 

I3 

I3 I4 ll I2 

I4 

I4 I3 I2 ll 


2 X li = 2 Vi, [2 X 2] = Ip + I 2 + I 4 , {2 X 2} = I 3 , 

where [n x n] denotes the symmetric and {n x n} the anti-symmetric part of the product 


n X n. 


The Clebsch-Gordan coefficients for the products 2 x are then with 


li^ 


51 

5 2 


2 and 


i = 1 i = 2 i = 3 


i = A 


( Sit \ f S2t \ f S2t \ 
\ S2t J \ Sit ) -Sit ) 


Sit 

-S2t 


~ 2 
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and for the product 2x2 for 


-Sl 

-S2 


h 

t2 




2 ; 


{siti + S2t2)/y/2 ^ Xl and {sit2 + S2ti) / \/2 ^ I2 , 
{sit2 - S2ti)/\/2 ~ I3 and (siti - S2t2)/V2 ~ I4 . 
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Figure 1 : Scatter plot of the correlation between sin^ 6^23 and |f/e3p for class B in case 
of normal hierarchy. All matrices associated with a zero entry in hi^ generate the same 
behavior. 



Figure 2 : Same as above for class B. All matrices associated with a zero entry in hi2 
generate the same behavior. 
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mfo/Amjoi 



Figure 3 : Scatter plot of \Uez\ against the ratio of the smallest neutrino mass squared 
divided by Am^ for class E, cf. Eqs. ( 1291311 ) . 
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